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ABSTRACT 

Visibility-visibility correlation has been proposed as a technique for the estimation 
of power spectrum, and used extensively for small field of view observations, where 
the effect of w — term is usually ignored. We consider power spectrum estimation 
from the large field of view observations, where the w — term can have a significant 
effect. Our investigation shows that a nonzero w manifests itself as a modification 
of the primary aperture function of the instrument. Using a gaussian primary beam, 
we show that the modified aperture is an oscillating function with a gaussian en- 
velope. We show that the two visibility correlation reproduces the power spectrum 
beyond a certain baseline given by the width, U w of the modified aperture. Further, 
for a given interferometer, the maximum U w remains independent of the frequencies 
of observation. This suggests that, the incorporation of large field of view in radio 
interferometric observation has a greater effect for larger observing wavelengths. 
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1 INTRODUCTION 

The directly observed quantity in radio interferometry is the complex visibility, measured as a 
function of baseline and frequency. Apart from being the building blocks for radio imaging (Per- 
ley et al., 1989), the observed visibility data can be used for the estimation of various statistical 
properties of the radio signal. Visibility based estimators have been widely used to estimate the 
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power spectrum of intensity fluctuations of the observed radio signal (Bharadwaj & Sethi, 2001; 
Bharadwaj & Ali, 2005; Dutta et al., 2009a) and have also been proposed to be a viable probe 
of Bispectrum (Saiyad Ali et al., 2006; Dutta et al., in preperation). For the Cosmic Microwave 
Background Radiation (CMBR) signal, the formalism Hobson & Maisinger (2002) shows the ef- 
fectiveness of using Maximum Likelihood Estimators (MLE) in the Visibility space for direct 
estimation of the angular power spectrum. A large number of radio telescopes, like the presently 
functioning GMRT 1 , upcoming MWA 2 and SKA 3 aim to probe the universe through the red- 
shifted 21 cm observations. The HI power spectrum is an unique observational probe of the dark 
ages (Bharadwaj & Sethi, 2001), the epoch of reionization (Datta et al., 2007) and the era, post- 
reionization (Bharadwaj et al., 2009; Saiyad Ali et al., 2006). The visibility based estimator has 
advantages over the image based methods for the estimation of power spectrum. Apart from the 
fact that it uses the data in its raw form, visibility based method is naturally more useful in cases 
of incomplete/sparse u-v coverage. The method has also been used on smaller galactic scales for 
estimation of HI power spectrum of a turbulent ISM (Dutta et al., 2008, 2009a,b) and continuum 
power spectrum of Supernova remnant (Roy et al., 2008). 

The definition of visibility V is often approximated as a 2D Fourier transform of the sky inten- 

FT 

sity distribution i.e, V(w, v) ^ 1(1, m), where (I, m) denote the two angular coordinates in the sky 
and (u, v), the corresponding Fourier conjugate variables. While, this is true for small field of view, 
and co-planar array distribution, the actual 3D visibility as is observed in radio interferometers, do 
not satisfy this simple relationship. For most practical purposes, at high observing frequencies, this 
approximation works reasonably well for the existing radio interferometers like the GMRT , VLA 
etc. However, upcoming radio telescopes aiming at wide field imaging will require inclusion of 
the full 3D effect. The same problem is also faced in imaging for large field of view observations 
or observations with non-coplanar baselines and is tackled by techniques of mosaicing (Perley et 
al., 1989; Frater & Docherty, 1980; BCornwell & Perley, 1992; Cornwell et al, 2008). 

However the problem of estimating power spectra using the measured three dimensional visi- 
bility V(u, v, w), has not been addressed. 

In this letter we investigate the effect of using the three dimensional visibility in the estimation 
of power spectrum and hence find a possible justification of the simpler two dimensional approach 
used extensively earlier. 
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http://www.gmrt.ncra.tifr.res.in/ 

http://www.mwatelescope.org/ 

http://www.skatelescope.org/ 
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2 EFFECT OF w — term ON THE APERTURE FUNCTION 

The direct observable in the radio interferometric observations is the complex visibility V 3D (U, v) . 

— * 

For a pair of antennae separated by d, with each antenna pointing along the direction of the unit 
vector k (referred to as the phase center) we have 

V 3D (U, v) = J dnne 2n0 < h - k) A(h - k, v)I(n - k, v) (1) 

where h denotes the unit vector to different directions of the sky, baseline U — d/X, A(n — k, v) 
denotes the primary beam and I(fi — k, v) is the specific intensity. Writing U = U± + wk, where 
U± is a 2D vector, and defining rl — k = 9,we have, for \9\ <^ 1, 9- k = 0, implying that ^is a 2D 
vector. In this limit 9 gives the position of any point on the sky with respect to the phase centre in 
a 2D tangent plane. This is known as the flat-sky approximation. The term wk quantifies deviation 
from this. 

In the 2D approximation we have 

V 2D ± , u) = J d 2 9 e 2 ™°^A(9, v) 1(9, v) (2) 

Writing the specific intensity 1(9, v) as 1(9, v) — I v + 51(9, v), where the first term is a constant 
background and the second term is a fluctuation, we have 

V 2D (U±, v) = I U A(U ± , v) + A(U ± , v) <g> SI(U ± , v) (3) 

where tilde represents a Fourier transform and ® denotes a convolution. 

The aperture function A(U±, v), peaks at U± — and has a finite width. Hence, we shall retain 
the second term in all subsequent discussions. 

We note that, ignoring the w — term leads to a simplification of the expression for visibility 
and in the 2D approximation, V 2D (U±, v) is the Fourier Transform of A(9, v)5I(9, v). Hence, we 
have, 

A(9, v) 51(9, v) = J d 2 U' ± V 2D (U' ± , v) e'^v'J (4) 
Substituting in equation (1) we obtain 

V 3D (U, u) = J d 2 U' L K(U, U' ± , v) V 2D (U' ± , v) (5) 
Where the kernel K(U, U' ± , v) is given as 
K(U,U' ± ,v) = J dn h e- 2 ™ (u' ± -u H n-k) 

= 4n 3o (2n\U' ± -U\), (6) 
with j denoting the Spherical Bessel function. 
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Figure 1. Modified aperture A(U± , w) plotted as a function of U± for two different values of w, (A) w = 400. A and (B) w = 12. kA at v = 150 
MHz. Solid line in B shows the gaussian envelope. 



Defining a quantity A(U, v) as 

A(U, v) = 4tt j d 2 U' ± Jo (2n\U' ± - U\)A(U' ± , v). (7) 
The 3D visibility takes the form 

V 3D (U, u) = j d 2 U' ± A(U - U' ± ) 5I(U' ± ) (8) 

It is to be noted that the 'w' dependence of V 3D (U, v) is translated to the function A(U — U±). 
This can be regarded as a modified aperture function. We investigate the nature of the modified 
aperture A(Uj_,w) as a function of U± at different values of w. We have assumed the primary 

-* r u 2 1 

aperture A(U±) to be a gaussian, exp — ^ja , of width U , and evaluated the integral in Eqn. (7) 
numerically. This indicates that A(U, v) has an implicit dependence on U . For an antenna of 
diameter D, U can be approximately written as U ~ D/X for observing wavelength A. We use 
Uq = D/X with D = 45 m (specifications of the GMRT) for the subsequent discussion. This 
corresponds to a field of view of 4° at 150 MHz. We shall discuss the effect of larger field of view 
later. 

Figure 1 shows the variation of the modified aperture A(U±, w) as a function of U± for two 
representative values of w (A: w = 400 A , B: w = 12 kA) at frequency v = 150MHz. 

For w = one reproduces the primary gaussian aperture trivially. We also note that for 
id <C (Vi the gaussian profile is still maintained (Figure 1:A). However for large values of w 
the aperture function manifest oscillations (Figure 1:B). The period of these oscillations is found 
to be sensitive to w (decreasing as w increases). The envelope of the modified aperture is also a 



gaussian, C exp 



u 2 

2U 2 



, with the parameter U w being a measure of it's dispersion. Ignoring the 



effect of oscillations in A, we note that the 3D formalism can be recast in the same form as it's 
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Figure 2. U w is plotted as a function of w for different central frequencies of GMRT. Horizontal solid line corresponds to the maximum possible 
baseline. 
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Figure 3. Modified aperture A(U± , w) plotted as a function of U± for w = 120. k A at v = 1420. MHz. 

2D counterpart with U w taking the role of U . Hobson & Maisinger (2002) have obtained a similar 
result assuming small field of view, where, they have shown that the effect of w-distortion can be 
considered as turning the primary beam into a complex gaussian. 

We next investigate the w and frequency dependence of U w . For a given frequency at large 
w, U w is found to increase linearly with w, i.e, U w ~ m(u)w. This implies that the w — term 
effectively broadens the aperture of the instrument. Figure 2 shows the variation of U w with w 
for different frequencies in log-log scale. We note that, the slope m(u) (as represented by the y- 
intercept in the Figure 2) determines the effect of w — term for increasing values of w. A small 
value of m(u) implies a slow increase of U w with w and the effect of the w — term is less. m{y) 
is found to fall off as ~ l/u with frequency. This indicates that the departure from the flat-sky 
approximation is more pronounced at the lower frequencies. Redshifted 21 cm line observed at 
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Hence, one may expect the w — term to have 



1420(MHz) _ j 



frequency v probes the redshift z = 
a greater effect while probing higher red- shifts. 

21 cm line has been used to study the ISM dynamics of the nearby galaxies (z ~ 0). Fig- 
ure 3 shows the modified aperture function for frequency 1420 MHz and w = 120 kA (this being 
the maximum U for GMRT like arrays). At this frequency, for GMRT, U = O.lkA, whereas 
U w \ w =u max — 0.18kA. It follows that the flat-sky approximation can be safely used if the largest 
length scale probed, corresponds to a U± 3> 3 U w \ w =u max - 

Till now, we have investigated the effect of w — term using D = 45 m, which corresponds to 
an field of view of 4° at 150 MHz. We estimated U w \ w =u max assuming D = 4mtoD = 45m 
at v = 150 MHz. At D = 4 m, (which corresponds to the largest proposed field of view, 45°, of 
SKA), U w | w =u max — 4 kA for U max — 120 kA. We have also observed that, U w \ w= u max 
This indicates that the effect of w — term is more pronounced for larger field of view, as expected. 



VISIBILITY CORRELATION AND POWER SPECTRUM ESTIMATION 

The power spectrum, P(U±, Au) of the intensity fluctuations 51 on the sky is defined as 

(8I(U±, vi) SI(U' ± , v*)) = P(U ± , Au) 5l(U ± - U' ± ) (9) 
where Av = \v± — We shall be considering, for simplicity, v x = 1^2 in all subsequent discus- 
sions and hence drop the Av dependence of the power spectrum. 
We define 

V 2 3D (U a , U b ) = (V 3D (U a ) V 3D *(U b )). (10) 

Using Eqn. (8) we obtain 

V 2 3D (U a ,U b ) = J d 2 Ul J d"u' ± A(U a - U'i) 

A*(U b -U' ± )(5I(Ul)5r(U' ± )) (11) 
Using the definition of the power spectrum and considering the correlation at the same base-line 

— * 

U, this simplifies to 

V? D (U) = J d 2 U' ± \A(U-U' ± )\ 2 P(U' ± ) (12) 

Noting that the effect of the w — term is contained in the modified aperture A we can retrieve the 
2D estimator V 2 2D used earlier Dutta et al. (2009a) by replacing A with A. Hence, 

V 2 2D (U ± ) = J d 2 U' ± \A(U ± - U' ± )\ 2 P(U' ± ). (13) 

We shall now discuss the effect of the w — term on the estimator defined in Eqn. (12). 
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Figure 4. V^ D as a function of U± for w = (dot-dash) and w = Umax (dash) at u = 150 MHz, assuming P(Uj_) = U±~ 2 . We also plot 
P(U±) = U±~ 2 (solid line) for reference. The vertical arrows show the U± value above which the power law is recovered at 1%. Note that the 
plots are given arbitrary offset for clarity. 



Figure 4 shows V% plotted as a function of U± for two values of w, (w — and w = U max ) 
at v — 150 MHz, assuming P(U±) = U±~ 2 . We have chosen U max = 12k\ (this being the largest 
baseline for the GMRT at 150 MHz). We have shown the power spectrum P(Uj_) = U±~ 2 for 
comparison. For large values of w, V^ show oscillations for U± < U w , which arises due to the 
oscillatory nature of \A(U)\ 2 . 

We find that V^ faithfully recovers the power law Uj_~ 2 (at 1%) for U± greater than a certain 
value. This value is found to be 3 x (y/2U ) for w = and 3 x (y/2U w ) at w — 12 kA. Hence, a 
non-zero w — term changes the U± value beyond which the power spectrum estimation would be 
valid. 

The quantity of interest in power spectrum estimation using the the radio interferometric ob- 
servations used earlier (Bharadwaj & Ali, 2005; Dutta et al., 2009a) is 



where, p(w) is a normalized probability distribution of w. The function p(w) is specific to an 
observation as well as to the array configuration of the interferometer. Hence, it is difficult to 
make a general quantitative statement regarding the effect of w — term in £. Since, for a given 
w, the largest baseline above which V~2 D (U± : w) ~ P(U±) is U w , we can qualitatively state that 
£ gives a good estimation of the power spectrum for U > 3 U w \ w =u max - ^ i s important to note 
that, for a specific array configuration, U w \ w =u max l% independent of the frequency v, whereas 
Umax oc v (Table 1). Hence, the U± range amenable for power spectrum estimation is larger for 
large observing frequencies. 




(14) 
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Table 1. The effect of w — term, comparison between various quantities. 
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0.18 









Table 2. Relevant U± (k A) values at different frequencies. 

3 DISCUSSION AND CONCLUSION 

We have studied the effectiveness of the widely used power spectrum estimator in presence of a non 
zero w. The w — term is found to affect the visibility correlation estimator through a modification 
of the aperture function (Table 1). The effect is more pronounced for the lower frequencies, where, 
for a given interferometer, the baseline range over which V 2 3D reproduces the power spectrum, 
is reduced. This restricts the largest possible length scales that can be probed using the 21 cm 
radiation from the epoch of reionization (20 < z < 6). However, for the observations of the 
nearby universe (i.e, v = 1420 MHz), w — term effect is not too significant (Table 2). 

The two-visibility correlation defined in Eqn. (12) has a positive noise bias associated with it, 
which can even exceed the signal in radio interferometric observations. To reduce the effect of the 
noise bias, one may follow the method used in the 2D analysis (Begum et al., 2006; Dutta et al., 
2008, 2009a) considering correlation of the visibilities at two nearby baselines (Eqn. (11), with 
Ub = U a + AC/). For simplicity, we have considered visibility correlations at the same baseline 
(Eqn. (12)). However, we note that an alternative method for reducing the noise bias is to observe 
the same field repeatedly for the same baseline configuration. 

Redshifted 21 cm observations allow us to probe the HI distribution of the universe over con- 
tinuously varying redshifts by tuning the frequency of the radio observations. By correlating vis- 
ibilities at different frequencies it is possible to do a tomographic study of the HI distribution, 
thereby probing the 3D power spectrum. In this work we have restricted ourselves to visibility 
correlations at the same frequency of observation. Future investigations may consider the effect of 
a nonzero 5v on power spectrum estimation. 
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